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We propose a self-consistent generalized quantum master equation (GQME) to describe electron 
transport through molecular junctions. In a previous study [M.Esposito and M.Galperin. Phys. 
Rev. B 79, 205303 (2009)], we derived a time- nonlocal GQME to cure the lack of broadening effects 
in Redfield theory. To do so, the free evolution used in the Born-Markov approximation to close the 
Redfield equation was replaced by a standard Redfield evolution. In the present paper, we propose 
a backward Redfield evolution leading to a time-local GQME which allows for a self-consistent 
procedure of the GQME generator. This approach is approximate but properly reproduces the 
nonequilibrium steady state density matrix and the currents of an exactly solvable model. The 
approach is less accurate for higher moments such as the noise. 



I. INTRODUCTION 

Recent advances in the experimental capabilities for 
constructing molecular junctions and measuring their re- 
sponse to external perturbations create new challenges 
for an adequate theoretical description of open quantum 
systems far from equilibrium.—^— In molecular junctions, 
most of the interesting applications are concerned with 
either strongly correlated or resonant tunneling regimes 
(or both). In these situations, conventional perturbation 
theory and/or effective mean-field theories (e.g. Hartree- 
Fock, GW, DFT) are inapplicable and may lead to un- 
physical predictions^^ 

In molecular junctions, molecules are sensible to pro- 
cesses such as oxidation, reduction, and excitations. This 
makes the description of transport in the language of 
the many-body states of the isolated molecules a con- 
venient tool for molecules weakly coupled to the con- 
tacts. Furthermore, sophisticated quantum chemistry 
methods use many-body states to describe molecular 
electronic structurcsj^ The dressed state representationfi 
often employed in quantum chemistry, is another ex- 
ample of many-body states formulation. The standard 
nonequilibrium Green function (NEGF) approach^^— is 
formulated in the language of elementary excitations and 
is therefore not well suited for a many-body states de- 
scription. A suitable alternative, the Hubbard NEGF 
approachfiir— still displays inconsistencies for low level 
of approximation that need to be resolved (see Ref. [l^ 
for a discussion). 

The simplest approach to transport at the molecu- 
lar state level is the Redfield quantum master equation 
(QME) originally developed in the context of NMRM 
and later in many other fields^ '^^'^^ including transport 
through quantum junctionsJ^^— It is derived using the 
Born-Markov approximation in order to get a closed evo- 
lution equation for the reduced density matrix. The ro- 
tating wave (or secular) approximation (RWA) ^^i^^ is of- 
ten invoked when the molecular levels are well separated 
so that the molecule Bohr frequencies evolve fast com- 



pared to the relaxation time scale induced by the cou- 
pling to the contact. This leads to a QME which has a 
simple physical interpretation in the molecular eigenba- 
sis. Indeed, populations satisfy a rate equation with rates 
given by Fermi golden rule (Pauli equation) and the co- 
herences are each independently exponentially damped. 
As a result, in nonequilibrium steady-states, all coher- 
ence effects (in the molecular eigenbasis) are lost. In 
its general form, the Markovian Redfield equation pre- 
dicts steady-state coherences that can be inaccurate and 
even sometimes lead to unphysical resultsJ^i^"— Fur- 
thermore, broadening effects are totally absent from the 
Redfield descriptioui^"— . In order to overcome some 
of these difficulties, we proposed in our previous work2£ 
a generalized QME capable of predicting the broaden- 
ing of the molecular levels in an approximate way. The 
key idea was to replace the assumption of free molecular 
evolution inside the kernel (which is second order in the 
molecule-contact interaction) used in the Born-Markov 
approximation by the Redfield evolution. The resulting 
equation was however non-local in time. Here we extend 
our consideration by proposing a kind of time-reversed 
Redfield evolution that leads to a time local generalized 
QME. This procedure is rewarding because it allows to 
formulate a practical self-consistent scheme to calculate 
the generator of the QME. 

The structure of the paper is the following. Section 
im introduces a model of molecular junction. In section 
lllli we briefly review our generalized QME derived ear- 
lier and introduce the self-consistent scheme. Section ITVl 
presents numerical examples and compares the results 
obtained within the scheme to a numerically exact ap- 
proach. Conclusions arc drawn in section IVl 



II. MODEL 

We consider a molecule (M) coupled to two metal con- 
tacts (L and R) each at its own equilibrium. The Hamil- 
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tonian of the system is 



(1) 



K=L,R 



The contacts are assumed to be reservoirs of free charge 
carriers 



(2) 



keK 



where K = L,R and (cfc) are creation (annihilation) 
operators for an electron in state k. 

The molecular Hamiltonian is represented in terms of 
the many-body states {|M)} of the isolated molecule 

Mi,M2 Mi,M2 

(3) 

where Xmi,M2 projection (or Hubbard) operators. In 
particular, in the eigenbasis representation, the molecular 
Hamiltonian reads 



Hm — ^ EmXm,m 

M 



(4) 



where Em are the eigenenergics. Note that the molecular 
many-body states \M) are characterized by all the rele- 
vant quantum numbers describing the state of an isolated 
molecule or a dressed state. 

The coupling between the molecule and the contacts is 
introduced in the usual way with hopping terms for the 
electrons moving from the contact to the molecule and 
vice versa 

^KM =Y.Y. {yMkX^M^k + VumcIXm) (5) 



keK M 



Here 



(6) 



is a transition from the molecular state Mi = Ni Si to the 
molecular state Mf = Nf s / . The number of electrons 
on the molecule is Nf = Ni — 1 and Si (s/) is the set of all 
the quantum numbers characterizing the molecular state 



in the charging block Ni {Nf), so that Xm = X 



M},Mi 



III. GENERALIZED QUANTUM MASTER 
EQUATION 

In this section we start by reviewing the derivation of 
the generalized QME introduced in our previous work. 
We then introduce our new way of closing the master 
equation by using a time-reversed effective evolution, and 
demonstrate how the resulting equation can lead to a 
self-consistent scheme to calculate the propagator of the 
molecular density matrix. 



Exact equation of motion 

We start by deriving exact equation of motion (EOM) 
for the quantity 

{X2i{t)) = Tr [e+'^(*"*«)X2ie-^^(*-*°)/5o] = (Ti2(i) 

(7) 

where |1) and |2) are many-body molecular states, <o is a 
starting point of the evolution usually taken at the infi- 
nite past to —oo, and po is the initial density operator 
for the whole system (molecule and contacts). We note 
that 



ai2it) = {l\TTK[m] |2) 



(8) 



is a matrix element of the reduced density operator of the 
molecule a at time t which is obtained by tracing out the 
contact degrees of freedom from the full density matrix. 

The Heisenberg EOM for the Hubbard operator in 
Eq.© yields an exact EOM for the reduced density ma- 
trix (see Ref. [2^ for details of derivation) which in molec- 
ular eigenbasis reads 



i?2)^12(t)-EE / 

M s •' - 



dt' 



+ ('^W2-lsa(0-^Jvi(*'))^A1(Af2-ls,2)(*' ^ ^) 
~ ^(2,7V2 + 1 s)M ^ i'){^M [t') ^l,N2 + l s (0) 
+ (^]/l(*')-^W2-ls,l(*)>^A1(Ar2-l.s,2)(*' ^ 

- (X2,Ari + ls(<)X]^(i'))S^(ijVi + is)(*' - 



-(-1) 



Mil,Ni + 



Here, M. is defined in Eq.®, S^'^^^ are the greater 
(lesser) projections of a self-energy due to the coupling 
to the contacts 

^MiM2{ti,t2) = Y VMik gk (n , T2) VkM2 (10) 

K=L.,RkeK 

where gk{Ti,T2) = — z(TcCfc(ri) c|(r2)) is the free elec- 
tron Green function (T^ is the contour ordering opera- 
tor). Dissipation matrix F is introduced in a usual way 

Tm,M2 {E) - I [^'m^M2 (E) - ^M,M2 (E)] (11) 

and within wide band approximation assumed below is 
energy independent. 

Eq. ^ is exact but does not have a closed form in terms 
of the reduced density matrix ai2. Indeed, its right hand 



(9) 
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side contains two-time correlation functions of Hubbard 
operators. Attempt to write EOM for the latter would 
lead to expressions having in their right hand side corre- 
lation functions of higher order involving both, molecular 
Hubbard operators and creation (annihilation) operators 
for the electrons in the contacts, at different times. This 
would generate an infinite hierarchy of EOMs. 



Generalized time-nonlocal QME 

The usual strategy involves closing the system of equa- 
tions by approximately representing high order correla- 
tion function in terms of correlation function(s) of lower 
order. In particular, closing Eq.® for the case when M 
are effective single particle orbitals at the first step (i.e. 
reducing the two-time correlation function in the right 
hand side to a single time average) is known as the gener- 
alized KadanofF-Baym ansatz (GKBA).— In our previous 
work^, we introduced the analogue of the GKBA in Li- 
ouville space where M are molecular many-body states: 

* E [SkM,M,it-t'){XM,MAt')Xl{t')) (12) 
Ml, Ah 

A similar expression holds for {xl^{t')Xi2{t)). Here 
gr(a) g^j.g j-g^aj-ficfj (advanced) Green functions in Liou- 
ville space 

^[2,34(0 = -ie{t) < A2i/B|e-*^*|X43/5^' > 

= -ie{t) < X2l|e-*^^"*|X43 > (13) 
^?2,34W = -lO{-t) « X^JbW^'\Xi2P'^ > 

= ie{~t) < XM\e'^'t'f'\Xi2 > (14) 

where = jf^ff^ is the equilibrium density operator 
for the contacts and Ib is unity operator in the contacts 
subspace. Using (fT2|) in ([9|), leads to a time- nonlocal 
generalized QME given by Eq.(35) in Ref. [H. 



Generalized time-local QME 

We now propose an alternative way to close Eq.® by 
reducing the two-time correlation function in the right 
hand side of ^ to the later rather than earlier time. In- 
deed, since the whole system evolution is time-reversible 
(one has to be careful with spin and magnetic field, 
though), and since the reduction of a two-time correlation 
function to a single-time average is an approximation, the 
reduction of the correlation function to a later time does 
not a priori seem to be worse than to the earlier one. 



The two-time correlation function can be exactly ex- 
pressed in Liouville space as 

{X,2{t)XUt'))^ 

9{t - t') <C 134/5 |e-'^(*'-*) \p{t)Xi2 » (15) 

+ e{t' - < xljB\e-'^^'-''^xl^p{t') » 

Then using the projection superoperator introduced in 
Ref.Si 

l^Mi.AhPB ^-^XmimJbI (16) 

Ml, Ah 

one gets the alternative ansatz 
{X,2it)xUt')) « 

* E [Sl2,AhAhit-t'){XMiAh{t')Xl{t')) (17) 
Mi,M2 

^{X,2it)Xl^MM^M.M2,3iit-t')_ 

and a similar expression for {xl^{t')Xi2{t)). Here we 
introduced the retarded and advanced Green functions 
on the Keldysh anti-contour— 

^^2,34(0 = -iO{-t) « X2liB\e-'^'\X43p'^ » 

= -~ie{-t) < X2i\e~' "^'"^IXas > (18) 

^?2.34(i) = < l34/B|e'^*|ll2PB' > 

= te{t) < X34|e* '^^"'\Xi2 > (19) 

where ^Ceff is the effective Liouvillian generating the 
time-reversed evolution. Its connection to the effective 
Liouvillian C^f / generating forward-time evolution is22ii2S 

^£ = C* (20) 

Note difference in sign with Refs. [28lj29l which is due to 
our definition of evolution operator as e~' ^"tf* instead 
of e'^"//* in Refs.Hilli. 

Applying (flT)) to ([9]) leads to a time-local version of 
the generalized QME 



_ai2(t) =E'^^12,34'T34(i) (21) 
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where in the molecular eigenbasis 



dt 



-^12, 34 -iiEl - E2) (5l,3 $2,4 - / 

S{2A)Mi'^) ^A1(l,3)(~*) ~ ^f2,4)x(*) ^!m(1,3)(~0 

+ <5l,3 (^Ar^.Ari X! [^°Af2-l s,4)Al(0 ^A<(iV2-l s,2) (^*) 
s 

+ ^(2.^2 + 1 s)A4 ^A4(4,JV2 + 1 s) (^0 
- ^2.4 (5^1,^3 X! [^("3,^1 + 1 s)7w(*) ^A^(l,iVi + ls)(~*) 

+S^i_ls,l)Al(0 ^!M(Ari-ls,3)(~0 } (22) 



Note that this expression for the effective Liouvillian in 
an arbitrary basis will differ from Eq. only in the free 
evolution term. 



If the QME, Eqs. ^ and ([22]), is used with the free 
molecular evolution instead of the effective dynamics in 
Eqs. ()18p and (fT9|) . one gets the standard Markovian Red- 
field QME. The RWA is justified when r,„„ < \ea - et\ 
{m,n S {a, &}) hi the molecular eigenbasis and consist 
in neglecting the nondiagonal of Tmn (here the standard 
effective single orbital formulation is used)^^. As a re- 
sult, coherences become decoupled from population in 
the molecular eigenbasis. The former die off at steady 
state and later obey a Pauli rate equation. When the 
RWA is not justified, coherence effects in the molecular 
eigenbasis cannot be neglected. 

It should be clear from our derivation that time-local 
form Eq. (PT|) of the generalized QME should have a sim- 
ilar degree of accuracy as the time-nonlocal version pro- 
posed earlier.— However, the present form is more suit- 
able for realistic calculations, but most important it nat- 
urally suggests a self-consistent procedure to evaluate the 
effective Liouvillian. Indeed, the effective Liouvillian, 
Eg.ip^. depends on the Liouville space Green functions, 
Eqs. ([T5| and p^ . which in turn depend on the Liouvil- 
lian through the connection (|20p . 



Expression for current 

The general expression for the time-dependent current 
at the contact-molecule interface K can be calculated 
within the non-equilibrium Hubbard Green function ap- 



proach a o-^^'^° 

lK{t) = 2Im /" dt' 

M,M' 

{XM{t)xi,,{t'))j:<,f^{t' ~t) (23) 
Using our ansatz, Eq. pT)) . leads to the expression 

/ + 00 
dt' Y (24) 

where M = {Nsf,N + Isi). 

Full counting statistics 



The theory of full counting statistics (FCS) was ini- 
tially proposed by Levitov and Lesovikj^i^^ and became 
popular in the molecular electronics community when 
shot noise in molecular junctions became experimentally 
measurable!^ The theoretical formalism for FCS within 
the QME was developed in Ref. Hi. Within the FCS the 
evolution operator is dressed by counting field(s) A which 
track the exchange of electrons between the molecule and 
the contact (s) 

keK M 

(25) 

Below we assume that counting starts at time to- 

The dressed Liouville equation for the total density 
matrix takes the form 



\px{t) cyi^[~iCx{t - to)] \p{h) > 



(26) 



As usual initial condition is assumed to be a direct prod- 
uct of the molecular and bath density matrices 

|p(to)»=|<T(to)/5B'» (27) 

so that dressed evolution of the system DM becomes 

\ax{t) exp[-i/:e//.A(t - to)] |a(to) > (28) 

The FCS is given by the generating function 

G(t,A) =«/|/)A(t) » (29) 

The long time limit of the logarithm of the generating 
function, 



5(A) = lim i In G(t, A), 

t-i-oo t 



(30) 
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provides information on the steady-state cumulants of 
the FCS 

In particular, the steady-state current is given by the first 
cumulant and the second cumulant yields to the zero- 
frequency shot noise. Using the spectral decomposition 
of the effective Liouvillian 

^e//,A = » >'lW < ^7(^)1 (32) 

7 

and since in the long time limit only one eigenmode vq (A) 
survives, we get 

^=-^^o(A) (33) 

S = ^^M>^) (34) 

We note that both ((24|) and ([33| provide the same steady- 
state current. 

IV. NUMERICAL EXAMPLES 

We now compare the results predicted by our new self- 
consistent approach with the standard Markovian Red- 
field equation and the nonequilibrium Green function re- 
sults. The simplest model which can be treated by these 
three methods, while providing information on both pop- 
ulations and coherences, is a non-interacting two-level 
bridge (TLB) between metallic contacts. 

The bridge part of the model has 4 many-body states: 
|0) = |0,0), \a) = |1,0), \b) = |0,1), and |2) = 
where \na, rib) indicates number of electrons Ua^b = {0, 1} 
on the level a and 6, respectively. The relevant single elec- 
tron transitions are: (0, a), (0,5), (&, 2), and (a, 2). 
They arc connected to the second quantized (single- 
particle) excitation operators used in standard GF ap- 
proaches by 

da = Xoa + Xb2 (35) 
4 = Xob - Xa2 (36) 

where {da,b) £^re the creation (annihilation) operators 
for an electron on level a and b, respectively. 

The molecular Hamiltonian, Eq.Q, in this case takes 
the form 

Hm ~ SaXaa + SbXbb + t (^Xab + Xba^ (37) 

The model was treated extensively within the QME 
approach 3^"— 

The TLB model is easily exactly solved using NEGF. 
The connection, Eqs. (j35|) and (|36p . provides partial in- 
formation about the many-body state populations but 



full information about coherences, which makes the com- 
parison between the different approaches meaningful. In 
particular, 

-tG<^{t,t)=aoo{t) + <7bb{t) (38) 
-iG<it,t)^aooit) + <Jaa{t) (39) 
^iG<b{t,t) aabit) (40) 

Inter-level coherence in the TLB is either due to hop- 
ing matrix element t or due to the coupling of the two 
levels via a common bath. The latter interference enters 
through the non-diagonal elements of the bridge-contact 
coupling matrix 

Ff, ^ VakVkbSiE - Ek) 

= Y.^ioa)kVkmSiE-ek) (41) 

keK 

= -Y^mkVkia2)S{E-ek) (42) 

where K = L,R. The matrix is energy independent in 
the wide-band limit assumed below. We intentionally 
perform simulations at low temperature T = 10 K to 
avoid broadening of the bridge levels due to artificially 
high values of the temperature. We focus on the relevant 
temperature regime for molecular junctions: F ^ ksT. 

We start our consideration from a simple far off- 
resonant case, where Tmn ^ \£a — £b\ {m,n G {a, 6}). 
Here our consideration is done in a local molecular ba- 
sis (t 7^ 0). The parameters of the calculation are the 
following: level positions Sa = 0.2 eV and Sb = 0.5 eV, 
inter- level hopping t = 0.1, escape rates for both levels 
into both contacts Ff„ = F^ = 0.1 eV (X = L, R), and 
level mixing due to coupling to the contacts F^ = F^ = 

0. 05 cV. The Fermi energy in the absence of bias is taken 
as zero, Ep = 0, and the voltage division factor is 1, 

1. e. fiL = Ef + \e\V and hr = Ep- The NEGF calcula- 
tion is performed on an energy grid spanning the range 
from -10 to 10 cV with step lO"* eV. Both the Rcdfield 
QME and our scheme are expected to work properly in 
this region of parameters. Figure [T] compares the re- 
sults of the Rcdfield QME and our new scheme to the 
exact NEGF results for the model. The main graph and 
the inset in FigH^ display populations, Eqs. ([55)1 and 
([55)1 . The real part of the coherence, Ea. (|l(I)) . is shown 
in FiglTjD. The imaginary part of the coherence is zero in 
this case (not shown). As discussed earlier;^ the Red- 
field QME approach misses the information about the 
broadening of the molecular levels. Our approach accu- 
rately recovers this information. Both populations and 
coherence are well reproduced by our GQME as well as 
the current-voltage characteristics (see Fig[T}:). Qualita- 
tively, the predictions of the Redfield QME approach are 
also correct here. 

We next consider the case of strong mixing due to 
coupling to the contact (s). Parameters of the calcu- 
lation arc t = 0, F^„ = 0.3 eV, F^„ = 0.048 eV, 
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FIG. 1: (Color online) Two-level bridge by Redfield QME 
(dash-dotted line, green), NEGF (dashed line, red), and self- 
consistent GQME (solid line, red) approaches. Shown are (a) 
probabilities, (b) coherence, and (c) current vs. bias. Far 
off-resonant local basis treatment. See text for parameters. 



FIG. 2: (Color online) Two-level bridge by Redfield QME 
(dash-dotted line, green), NEGF (dashed line, red), and self- 
consistent GQME (solid line, red) approaches. Shown are (a) 
probabilities, (b) coherence, and (c) current vs. bias. Case of 
strong mixing due to coupling to the contact (s). See text for 
parameters. 



rf, = 0.2 eV, rg = 0.1125 eV, T^, = Vg = 0.12 eV, and 
r^f, = = 0.15 eV. The other parameters are the same 
as in FiglT] Figure[2]comparcs the Redfield QME and our 
self-consistent GQME scheme with the exact NEGF re- 
sults. Our scheme still reproduces populations (Figl5^), 
coherences (Fig[2}D), and current (Figl2]:) accurately. The 
Redfield QME however fails to reproduce, even qualita- 
tively, the populations (see Figl2^) and the real part of 
the coherence (see Figl^b). This is due to coherence ef- 
fects through the bath which are not properly captured 
by the Redfield QME. 

Interference effects in molecular systems were ob- 



served experimentally for electron transfer— and molec- 
ular junction currents^ involving derivatives of benzene 
connected in meta or para position. They were also ex- 
tensively discussed in theoretical literature42ri^ Figure [3] 
shows current (a) and conductance (b) vs. bias, for a 
model system in which destructive interferences can be 
experimentally observed and which has been discussed 
earlier in the literature4^ This model is a two-level sys- 
tem with only one of the levels, a, attached to both L 
and R contacts. The other level, b, is coupled to the 
level a through a hopping element t. Level b is not di- 
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FIG. 3: (Color online) Two-level bridge by Redfield QME 
(dash-dotted line, green), NEGF (dashed line, red), and self- 
consistent GQME (solid line, red) approaches. Shown are (a) 
current and (b) conductance vs. bias. Resonant, Sa ~ £b, 
consideration in a local basis with only level a coupled to 
contacts, rf„ / and = F^^ = = (A' = L, R). See 
text for parameters. 



rectly attached to any contacts. As a result, the tunnel- 
ing electron has two possible pathways to be transferred 
from contact L to contact R via the system: one di- 
rectly through level a and the other by exploring level b 
on its way. Interference between the two paths leads to 
destructive interference in the transport characteristics, 
which reveals itself as a dip in the conductance. Parame- 
ters of the calculation arc Sa ~ £b = 0.25 eV, t = 0.05 cV, 
= ^1 = 0.2 eV and Tf, = Tf, = = 0. Figll. 
shows that destructive interference in conductance is ac- 
curately captured by our self-consistent GQME approach 
but not by the Redfield QME. 

Finally, we present the results of calculation for the 
zero-frequency shot noise. The shot noise was calculated 
using the PCS approach, which for the Redfield QME 
and our self-consistent scheme yields to expression ([M]). 
The PCS within NEGP was discussed in our previous 
publication;^ Note that the PCS for a two-level bridge 
within QME was also discussed in Ref. [s^ We intention- 
ally consider two independent levels (no coupling either 
within the system or through the bath) to demonstrate 



FIG. 4: (Color online) Two-level bridge by Redfield QME 
(dash-dotted line, green), NEGF (dashed line, red), and self- 
consistent GQME (solid line, red) approaches. Shown are (a) 
shot noise and (b) differential shot noise vs. bias. See text 
for parameters. 



the general problem that usual QME schemes have in 
representing higher moments. We assume that the level 
Ea = 0.5 eV is coupled symmetrically to both contacts 
rL = ^aa = 0.1 eV. The other level = 1.5 eV is 
coupled asymmetrically P^^ = lOP^ = 0.1 eV. Other 
parameters arc t = 0, Pf^, = P^^ = 0, and T = 10 K. 

Figure 2] compares the results of the Redfield QME 
and the self-consistent GQME with the exact results pro- 
vided by NEGF. It is known^ that the differential shot 
noise yields a two-peak structure for a symmetrically cou- 
pled molecule, with only a single peak observed for highly 
asymmetric coupling. One sees that the latter peak (cen- 
tered around 1.5 eV in FigHJa) is reproduced quite well by 
our self-consistent approach. The double well structure 
is however missed by both the Redfield QME and our 
approach (see peak centered aroundi 0.5 eV). Cumulants 
beyond first order can be very sensible to approximations 
and our level of theory, which neglects molecule-contacts 
correlations beyond second order, has to be improved to 
properly reproduce them. 

Summarizing, within our simple two-level model of 
molecular junction, the comparison of the proposed self- 
consistent GQME and the Markovian Redfield QME to 
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the exact results provided by the NEGF shows that the 
self-consistent GQME is highly accurate even when the 
Redfield QME fails qualitatively. We find that the self- 
consistent GQME will only fail in the region of parame- 
ters where the inter-level coherence due to the coupling 
to contacts in the molecular eigenbasis is bigger than 
the distance between the molecular eigenstates and is 
of the order of the molecule-contact coupling strength 
(i.e. exactly at resonance). As a result, in most relevant 
practical calculations of transport in molecular junctions, 
we expect our approach to be accurate. The ability of 
the proposed scheme to treat molecular transport in the 
language of many-body states of the isolated molecule 
and its ability to properly account for interference effects 
makes it a valuable practical tool for ab initio calcula- 
tions of transport in molecular junctions. 

V. CONCLUSION 

We present a practical scheme for molecular transport 
calculation. The scheme is based on a time-local gener- 
alized quantum master equation obtained by closing the 
exact EOM for Hubbard operators by employing a time- 
reversed evolution ansatz on the Keldysh anti-contour, 
similar to generalized Kadanoff-Baym ansatz introduced 
in our previous publication. We note that the approxi- 
mations involved in derivation of the time-local equation 
are essentially the same as for the earlier (more tradi- 
tional) time-nonlocal GQME. The time-locality of the 
GQME allows us to formulate a feasible self-consistent 
scheme to calculate the time dependent molecular density 
matrix and the current in term of molecular many-body 
states. We find that the convergence of the self consistent 
method for a simple two-level bridge model is achieved 



within two iterative steps. The results of the calcula- 
tion for TLB within self-consistent GQME are compared 
to Redfield QME approach, and to exact NEGF results. 
We demonstrate that our scheme (contrary to the usual 
QME result) properly captures populations and coher- 
ences. In particular, destructive interference effects in 
the molecular devices previously discussed in the liter- 
ature can now be properly described in the many-body 
states language, which makes the scheme a valuable tool 
for practical ab initio calculations. Current-voltage char- 
acteristics are reproduced with high accuracy. We are 
able to calculate the molecular device characteristics in 
the experimentally relevant regime of low temperatures. 
We find that the scheme breaks down in the region of 
the parameters where coherences in the system eigenba- 
sis (i.e. coherences introduced through non-diagonal el- 
ements of molecule-contact coupling matrix F) arc both 
bigger than the inter-level separation and are of the order 
of the escape rates (diagonal elements of the molecule- 
contact coupling matrix F). Thus, we expect the scheme 
to be a valuable tool for most practical transport calcu- 
lations in molecular junctions. 

The formulation of a scheme capable of reproducing 
higher moments of the full counting statistics in the 
many-body state language that goes beyond the Redfield 
QME approach^ is the goal of future research. 
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